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Ultraluminous Infrared Galaxies at 1.5 < 2; < 3 occupy dark 
matter haloes of mass ~ 6 x IO^^Mq 
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Abstract. We present measurements of the spatial clustering of ultralumi- 
nous infrared galaxies in two redshift intervals, 1.5 < z < 2.0 and 2 < z < 3. 
Both samples cluster strongly, with rp = 14.40 ± 1.99ft-^^Mpc for the 2 < 2 < 3 
sample, and tq = 9.40 ± 2.24/i~^Mpc for the 1.5 < z < 2.0 sample, making 
them among the most biased galaxies at these epochs. These clustering ampli- 
tudes are consistent with both populations residing in dark matter haloes with 
masses of ~ 6 x IO^^Mq. We infer that a minimum dark matter halo mass is 
an important factor for all forms of luminous, obscured activity in galaxies at 
z > 1. Adopting plausible models for the growth of DM haloes with redshift, 
then the haloes hosting the 2 < z < 3 sample will likely host the richest clusters 
of galaxies at z=0, whereas the haloes hosting the 1.5 < z < 2.0 sample will 
likely host poor to rich clusters at z—Q. 



1. Introduction 

When we look at the night sky, we see that galaxies seem to be arranged in a 
particular way. One might expect that galaxies would be distributed randomly, 
much as grains of sand would if you threw a handful across the floor, but instead, 
they seem to trace elegant structures; galaxy clusters are connected to each other 
by long filaments, interspersed with large voids, where few or no galaxies are 
seen. The drivers behind the formation of these 'large-scale structures' have 
been the subject of intense study and debate for over thirty years; how did the 
Universe go from being smooth and homogeneous just after the Big Bang to the 
clumpy, clustered Universe we see today? 

At the core of current theories for the formation of these structures is the 
premise that the evolution of the total mass distribution is described by the 
gravitational collapse of primordial density fluctuations, and that this evolu- 
tion is traced by the evolution of galaxies. Overdense regions, or 'haloes', are 
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predicted to undergo mergers to build haloes of increasing mass, with galaxies 
forming from the b aryonic matter in these haloes. This framework of 'biased' hi- 
erarchical buildup (jCole et alll200Q : Hatton et al 2003 ) has proven to be remark- 
ably successful in explaining several important aspects of galaxy and large-scale 
structure formation. 

This paradigm is however not without its problems. A good example of 
these problems concerns the evolution of massive (> 1O"M0) galaxies. We 
might expect that massive galaxies form slowly, with many halo mergers needed 
to build up the required larg e baryon reservoirs, and indeed some ga laxies do 
appear to form in this way (|van Dokkum et all Il999l : iBell et all l2004l ) . There 
i s however evidenc e that many massive galaxies may form at hig h redshift 
( Dunlop et all [199^ : iBlakeslee et all [20031 ) and on short timescales ( Ellis et al 
19971 ), direct ly counter to ear l y, 'naive' inodel predictions. Intriguingly, recent 



surveys (e.g. lEales et"alll2000l : Scott et al |2002| ) have uncovered a huge popula- 
tion of distant, IR bright sources that are plausible candidates for being rapidly 
forming, massive galaxies, however they are so numerous that even recent models 
have difficulty in explai ning their counts, and invoke a wide variety of solutions 
(e.g. iBaugh et all (|2005l 11. 

It seems likely therefore that there are strong, but subtle links between 
distant, IR/sub-mm bright galaxies, and the formation of large-scale struc- 
tures. Therefore, we need observations that relate the properties of these galax- 
ies with the underlying dark matter distribution. Motivated by this, we have 
used data from the S pitzer Wide Area Infrared Extragalactic Survey (SWIRE, 
Lonsdale et all ( 20031 )) to select large samples of distant Ultraluminous Infrared 
Galaxies (ULIRGs, Ljj. > IO^^Lq) and study their clustering evolution with 
redshift. We assume Hq = lOO/i km s"^ Mpc~\ O = 1, and Oa = -7. The 
results presented here were originally published in Farrah et al. ( 2006al lbl). 



2. Analysis and results 

To select high redshift ULIRGs, we use the 1.6/im emission feature, which arises 
due to photospheric emission from evolved stars. When this feature is red- 
shifted into one of the IRAC channels th en that channel exhibits a 'bump' 
( Simpson &: Eisenhardtlll999l : ISawickill2002l ). A complete discussion of the source 



selection and characterization methods is given in Lonsdale et al 2006 (in prepa- 
ration), which we summarize here. 

Our sources are taken from the ELAIS Nl and ELAIS N2 fields, and the 
Lockman Hole, covering 20.9 square degrees in total. We first selected those 
sources fainter than R=22 (Vega), and brighter than 400/iJy at 24/im. Within 
this set, we selected two samples that displayed a 'bump' in the 4.5//m and 5.8//m 
channels, i.e. where /s.e < /4.5 > /s.s > /s.o for one sample (the 'B2' sample) 
and where /a.g < /4.5 < /s.s > /s.o ^or the other sample (the 'B3' sample). This 
resulted in a total of 1689 B2 sources an d 1223 B3 sources. 

For both samples we used Hyper-z ( Bolzonella et alll2000l ) to estimate red- 



shifts, the results from which place most of the B2 sources within 1.5 < z < 2.0, 
and most of the B3 sources within 2.2 < z < 2.8. From the best fits we also 
derived IR luminosities and power sources; the requirement that the sources 
have /24 > 400;uJy demands an IR luminosity of > IO^'^Lq for all the sources. 
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with (most objects having) a starburst as the dominant power source, with star 
formation rates of > 2OOM0yr~^. Similarly, the presence of the 1.6/um feature 
demands a minimum mass of evolved stars of ~ lO^^M©. Both the B2 and B3 
sources are thus good candidates for being moderately massive galaxies harbor- 
ing an intense, obscured starburst, making them similar in nature to both local 
ULIRGs, and high-redshift SMGs. 

We me asured the angular c lustering of both samples using the methods 
described in Farrah et al. ( 2006al lbl). which we summarize here. We found that 
the levels of angular clustering seen in the three fields were consistent with 
each other to within 0.5a, and so combined the angular clustering measures for 
each sample over the three fields. To quantify the strength of clustering, we fit 
both datasets with a power law, u:{6) = Ai^6^~'^, where 7 = 1.8 and A^^ is the 
clustering amplitude; = 0.0125±0.0017 for the B3s and A^ = 0.0046±0.0011 
for the B2s. To convert these angular clustering amplitudes to spatial clustering 
amplitudes, we invert Limbers equation: 



^o(^) 



Hq j4(^ cC 



Ja dz 
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(1) 



where f{z) parametrizes the redshift evolution of tq, and: 
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W e derived dN/ dz from the photometric redshift distributions (jLonsdale et al 
20061 ). For the B2 sources this is a Gaussian centered at 2: = 1.7 with a FWHM 
of 1.0, and for the B3 sources this is a Gaussian centered at z = 2.5 with a 
FWHM of 1.2. The resulting correlation lengths are ro = 9.4 it 2.24/i~^Mpc for 
the B2 sources, and tq = 14.4 it 1.99/i^^Mpc for the B3 sources. 



3. Discussion 

To place these clustering results in context, we consider two models. The first 
parametrizes the spatial correlation function, ^, as a single power law in comov- 
ing coordinates, where the comoving correlation length, rQ{z), is: 

ro(z) = ro/(z),/(z) = (l + z)T-(3+^) (3) 
Here the choice of e determines the redshift evolution ( Phillipps et al. 19781 : 



Overzier et al|[2003l ). Several cases are usually quoted. First is 'comoving clus- 



tering', where haloes expand with the Universe, and e = 7 — 3; in this case 
clustering remains constant. Second is the family of models for which e > 0, 
for which clustering increases with time. Examples of this family include (a) 
'stable' clustering, for which e ~ (in this case the haloes are frozen in proper 
coordinates, (b) the predicted evolu tion of clustering of the overall dark mat- 
ter distribution, wh ere e ~ 7 — 1 (jCarlberg et al|[2000l ). and ro ^ 5 at z=0 



( Jenkins et a]|[l99^ ). (c) 'linear' clustering, where e = 1.0. A cautionary note 



to this is that detailed interpretations of clustering evolution from these models 
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suffer from several theoretical flaws ( Moscardini et allll998l : ISmith et ani2003l ). 
and so should be thought of as qualitative indicators rather than quantitative 
predictions. We therefore simply use the 'stable' and 'linear' models as indicators 
of the possible range of halo clustering amplitude with redshift. 

The second class of model comprises those in which comoving correlation 
lengths increase with increasing redshift. These models introduce 'bias', b{z), 
between the galaxies and the un derlying dark inatter . An example of suc h 
models are the 'fixed mass' models ( Matarrese et aCl997l : iMoscardini et allll998l ). 
which predict the clustering strength of haloes of a specified mass at any given 
redshift. In Figure [1] we plot the e model for dark matter, 'stable'and 'linear' 
epsilon models normalized to the B2 and B3 clustering strengths, the fixed halo 
mass models for halo masses of lO^^M©, lO-'^^M©, and IO^^Mq, the ro values 
for the B2 and B3 galaxies, and the spatial correlation lengths of other galaxy 
populations taken from the literature. 

With the the uncertainties described earlier firmly in mind, we use Figure 
[T] to explore the relationships between our samples, the underlying dark mat- 
ter, and other galaxies. Both the B2s and the B3s are strongly clustered, with 
correlation lengths much higher than that predicted for the overall DM distri- 
bution. Both B2s and B3s cluster significantly more strongly than optical QSOs 
at thei r respective epochs, a nd B3s cluster more strongly than SMGs. Based 
on the Matarrese et al ( 19971 ) models, then we derive approximate la halo mass 
ranges of lO^^.T ^^^^ ^ iqU.i ^j^g ggg^ ^qIS.s < lO^^.Q ^.j^g 

B2s. Interestingly, halo masses comparable to th ese were recently der i ved fo r an 
independent sample of high-redshift ULIRGs bv iMagliocchetti et al.l ( 20071 ). 

The most interesting comparison is however betwee n the two sarnples t hem- 
selves. The clustering evolution of QSOs with redshift ( Groom et al. 20051 ) may 
mean that there is a 'minimum' host halo mass for QSO activity, below which 
no QSO is seen, of ~ 5 x lO^^M©. The correlation lengths for the B2 and B3 
samples are consistent with the same conclusion but for a ~ 6 X IO^^Mq ' mini- 
mum' DM halo mass. Taken together, these results imply that a minimum halo 
mass is a 'threshold' factor for all forms of luminous activity in galaxies, both 
starbursts and AGN. It is also interesting to speculate on what the host haloes 
of B2 and B3 sources contain at lower and higher redshifts. We might expect 
that a halo hosting a B3 source could contain an optically bright LBG at z ~ 4, 
followed by a B3 at z y 2.5, possibly ac c ompa nied by other (near-IR selected) 
star forming systems ( Daddi et al. 2004 . 20051 ). before evolving to host a rich 
galaxy cluster at low redshifts. The occupants of a halo hosting a B2 galaxy 
would however probably be different. We would expect that such a halo could 
contain an SMG at z ~ 2.5, and optically fainter LBGs at 4 < z < 5 (though 
probably not LBGs at z ~ 3). At lower redshifts such a halo might host a 
radio-bright AGN and or ERO at z ~ 1, and a (poor to rich) cluster at z = 0. 
We conclude that ULIRGs at z > 1.5 as a class likely signpost stellar buildup 
in galaxies in clusters at z = 0, with higher redshift ULIRGs signposting stellar 
buildup in galaxies that will reside in more massive clusters at lower redshifts. 
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Figure 1. Comoviiig correlation l ength, rp, vs. redsh i ft. Other data are 
taken from Moscardini et al 1998; Ovcrzier et al 2003; Daddi et al.l 



2004; 



2005 



Blain et al2004; Ouchi e t al. 2004; A dclbcrgcr et al. 2005; Gr oom et al.l 
Georgakakis et al.ll2005l : lAllen et ai.l 2005. ^The 'Fixed mass' lines show the 



predicted clustering amplitude of haloes of a given mass at any particular 
redshift, whereas the e lines show the predicted clustering amplitude of an 
individual halo for three halo growth models, described in the text. The 'Sta- 
ble' and 'Linear' lines give a qualitative indicator of the range of how DM 
haloes may grow with redshift, and we have normalized 'Stable' and 'Linear' 
lines to the clustering amplitudes of the B2s and the B3's. The shaded regions 
therefore indicate what these haloes may host at lower and higher redshifts - 
the haloes hosting B3s may contain an optically bright LBG at z ~ 4 (upper 
green point), and grow to host very rich galaxy clusters at z=0, whereas the 
haloes hosting B2 sources may contain optically fainter LBGs at 4 < z < 5, 
SMGs at z ~ 2.5, radio-bright AGN (upper pink triangle) and (old) EROs at 
z ~ 1, and poor to rich clusters at z=0. 
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